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ABSTRACT 

Direct  heat  addition  of  flowing  gases  by  radiation  absorption  is 
considered.  In  the  IR  regime  the  coupling  between  a  laser  beam  and  a 
flowing  gas  is  being  modeled  numerically.  Implicit  time-iterative 
procedures  originally  developed  for  transonic  flows  are  being  adapted  to 
the  low  Mach  number,  low  Reynolds  number  regimes  of  interest. 
Two-dimensional  solutions  show  that  absorption  plasmas  can  exist  over  a 
wide  range  of  flow  speeds  and  that  buoyancy  is  a  dominant  factor  in  the 
forced  convection  flowfields.  In  the  microwave  regime,  an  experimental 
investigation  of  various  absorption  modes  is  underway.  An  overview  of  the 
experimental  design  and  initial  experimental  results  is  presented. 
Numerical  predictions  of  the  resonant  cavity  and  propagating  plasma  modes 
of  microwave  absorption  are  also  discussed. 


I. 


STATUS  OF  RESEARCH  EFFORT 


The  advanced  computational  algorithms  that  have  been  developed  for 
transonic,  external  aerodynamics  in  the  past  decade  are  being  adapted  for 
low  speed  flows  with  strong  heat  addition.  The  specific  area  of  focus  has 
been  on  the  computation  of  laser  gas-dynamic  interactions,  but  the 
techniques  are  equally  applicable  to  solar  energy  absorption  in  flowing 
gases,  arcjets  and  other  advanced  propulsion  concepts  as  well  as  to  the 
more  traditional  combustion  problems  that  occur  in  propulsion  environments. 
The  family  of  techniques  being  addressed  is  time  dependent  procedures. 

These  are  attractive  because  they  apply  to  either  viscous  or  inviscid  flows 
and  require  a  minimum  amount  of  artificial  viscosity.  In  addition,  time 
dependent  procedures  are  applicable  to  either  steady  or  unsteady  flows, 
making  them  candidates  for  studying  both  the  steady  flow  characteristics  of 
propulsion  environments  as  well  as  the  stability  characteristics  of  these 
flowf ields. 

The  extension  of  time  dependent  algorithms  from  external  transonic 
flows  to  internal  propulsion  environments  requires  advances  in  several 
directions.  First,  the  methods  must  be  modified  to  handle  very  low  Mach 
number  flows  (which  are  nonetheless  compressible  because  of  the  energy 
deposition)  as  opposed  to  the  transonic  flows  for  which  they  were  developed. 
Second,  they  must  be  extended  from  high  Reynolds  number  applications  to  low 
Reynolds  number  conditions  for  at  least  some  of  the  energy  deposition 
problems  of  interest.  The  laser  absorption  problem  in  particular  is  a 
phenomenon  that  takes  place  at  low  Reynolds  numbers.  Very  effective 
extensions  to  low  Mach  number  have  been  accomplished  as  discussed  below. 

Work  on  developing  an  effective  low  Reynolds  number  algorithm  is  still  in 
progress.  Finally,  the  consideration  of  strong  heat  addition  at  low  speeds 
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implies  that  the  effects  of  buoyancy  must  be  addressed,  and  substantial 
progress  has  been  made  along  this  line  also. 

The  equations  of  motion  for  a  flowing  gas  with  heat  addition  (by 
absorption  of  radiation,  I2R  heating,  or  combustion)  can  be  expressed  in 
vector  form  as, 


3Q  _  £E  +  3F  n  ,  3Ey  ,  3FV 
at  3x  3y  3x  ay 


(1) 


Here,  Q,  E,  F,  H,  Ev  and  Fy  are  vectors.  The  primary  dependent  variables 


are  contained  in  Q  while  the  remaining  five  vectors  are  (in  general) 
nonlinear  functions  of  Q.  These  vectors  are  defined  as, 


H 


and, 

E  m  R  ^,0^  ^  p  _c^Qiyr 
v  xx  3x  xy  3y 

p  ,  R  12a  +  R  12a 
v  yx  3x  yy  3y 

where  Qy  *®  (p,  u,  v,  T)^*,  and  Rjjy,  Ry^,  end  Ryy  ere  matrices 

containing  only  the  viscosity  and  the  thermal  conductivity.  For 
computational  purposes,  the  equations  are  transformed  to  generalized 
nonorthogonal  coordinates.  Convergence  to  a  steady  state  corresponds  to 
following  Eqn.  1  through  a  transient  from  some  arbitrary  initial  condition 
to  a  final  steady  condition.  We  begin  by  discussing  the  inviscid  terms  on 
the  left-hand  side  and  then  add  the  viscous  terms  on  the  right-hand  side. 


It  is  well  known  that  the  convergence  range  of  typical  time  dependent 
methods  becomes  slower  as  the  flow  Mach  number  is  reduced.  The  reason  is 
because  the  eigenvalues  of  the  system  become  increasingly  stiff  as  Mach 
number  is  reduced.  We  have  developed  two  methods  for  circumventing  this 
problem.  In  the  first,  the  time  derivatives  are  multiplied  by  some  matrix 
r  which  is  chosen  in  such  a  way  that  the  eigenvalues  remain  well 
conditioned  at  low  Mach  numbers.  This  matrix  preconditioning  procedure 
corresponds  to  modifying  the  physical  time  derivatives  so  that  the 
preconditioned  derivatives  provide  faster  convergence.  Our  first  attempts 
at  preconditioning  allowed  rapid  convergence  down  to  a  Mach  number  of  0.01, 
but  became  ineffective  below  that.  In  addition,  extension  to  real  gases 
with  variable  properties  showed  that  the  preconditioning  was  not  effective 
when  specific  heats  varied  strongly.  Later  refinements  of  the 
preconditioning  allowed  us  to  maintain  rapid  convergence  to  much  lower  Mach 
numbers  but  eventually  round-off  errors  were  encountered  in  the  pressure 
computations.  These  round-off  errors  effectively  prevent  the  use  of  the 
procedure  at  Mach  numbers  below  about  10"^.  Testing  of  the  refined 
preconditioning  scheme  for  real  gas  equations  of  state  has  not  yet  been 
done,  but  it  is  possible  that  preconditioning  will  now  be  effective  for 
real  gases.  These  results  have  been  summarized  in  Ref.  9. 

The  second  method  for  accomplishing  low  Mach  number  computations  has 
been  to  use  a  perturbation  expansion  of  the  equations  of  motion  to  obtain  a 
low  Mach  number  system.  In  this  low  Mach  number  system,  artificial 
acoustic  modes  were  introduced  that  provided  a  well-conditioned  eigenvalue 
system  at  any  Mach  number.  This  system  also  included  a  rescaling  of  the 
pressure  and,  thus  far,  all  inviscid  applications  have  shown  it  to  be  very 
effective.  We  have  obtained  successful  convergence  to  Mach  numbers  as  low 
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as  10~6  which  exceeds  the  low  speed  range  of  the  laser  propulsion  problem 
of  interest.  The  low  Mach  number  expansion  procedure  is  described  in  Refs. 
9  and  10. 

At  these  very  low  speeds,  buoyancy  becomes  a  significant  parameter. 
Buoyancy  introduces  a  source  term  into  the  momentum  equations  and  source 
terms  (as  opposed  to  sink  terms)  are  always  destabilizing  in  a  numerical 
solution.  The  effects  of  the  source  term  on  convergence  can  be  mitigated 
by  an  alternative  preconditioning  that  effectively  makes  the  equations 
become  stiff  again,  but  results  have  shown  that  adequate  convergence  can  be 
maintained  to  Mach  numbers  as  low  as  10-J4  in  argon  with  heating 
representative  of  that  in  laser  propulsion  environments.  Additional 
potential  methods  for  counteracting  the  effect  of  buoyancy  have  been 
developed,  but  not  implemented  yet. 

The  viscous  terms  also  have  a  different  effect  on  convergence  at  low 
Reynolds  numbers  than  they  do  at  high  Reynolds  numbers.  Implicit  schemes 
traditionally  use  approximate  methods  for  solving  the  large  sparse  matrices 
that  are  generated.  These  approximate  methods  introduce  errors  of  order 
At2  that  become  large  and  dominate  convergence  when  large  time  steps  are 
involved.  In  the  case  of  the  Euler  equations,  these  error  terms  represent 
products  of  lnviscid  derivatives  in  the  x-direction  with  inviscid 
derivatives  in  the  y-direction.  The  resulting  approximate  algor  ;hm  is 
unconditionally  stable.  Similarly,  when  only  the  viscous  terms  are  present 
the  error  terms  are  composed  of  viscous  derivatives  in  x  times  viscous 
derivatives  in  y.  Again,  the  approximate  algorithm  for  the  viscous 
equations  remains  unconditionally  stable.  Because  both  the  viscous  and 
inviscid  problems  are  stable,  It  is  logical  to  expect  that  combined 
viscous-inviscid  problems  would  also  remain  stable,  but  stability  analyses 
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and  numerical  experiments  demonstrate  this  is  not  the  case.  Approximate 
factorization  of  the  combined  viscous-inviscid  problem  also  generates 
cross-derivati v  error  terms  such  as  inviscid  derivatives  in  x  times 
viscous  derivatives  in  y.  These  additional  error  terms  are  destabilizing 
and  cause  the  algorithm  to  be  unstable  in  the  Reynolds  number  range  of 
interest.  Several  methods  for  circumventing  this  difficulty  have  been 
tested  though  none  has  yet  been  found  that  is  completely  successful. 

The  primary  direction  being  pursued  is  the  use  of  an  alternative 
formulation  of  the  energy  equation.  When  the  diffusive  terms  were  added  to 
the  form  of  the  energy  equation  that  was  used  for  inviscid  flow,  it  became, 


l£l  + 
at 


vn 

Yp°(3x  3y 


a_  ,  3T  +  a_  ,3T 

3x  3x  3y  3y 


(2) 


Here,  p^  is  the  perturbation  pressure  and  the  time  derivative  represents  an 
artificial  derivative  that  is  used  for  convergence  purposes.  The  quantity, 
pQ,  is  the  nominal  or  mean  pressure  in  the  flowfield.  Both  theoretical 
(stability  results)  and  experimental  (computer  solutions)  analyses  of  this 
form  showed  that  it  was  unstable  for  viscous  flows.  To  avoid  this 
instability,  the  energy  equation  was  modified  by  combining  it  with  the 
continuity  equation  and  changing  the  time  derivative  to  one  containing  the 
temperature  to  give, 


3T  + 
3t  + 


h  puT  *  h pvT 


3_  3T  +  3_  3T 
3x  3x  3y  3y 


(3) 


This  has  slightly  better  convergence  than  does  Eqn.  2,  but  it  still  remains 
unstable  at  Reynolds  numbers  of  about  100.  (Note  that  when  X  -  0,  both 
Eqns.  2  and  3  perform  very  well  in  inviscid  flows.) 
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In  full  vector  form,  the  perturbation  equations  corresponding  to  Eqn. 
3  can  again  be  written  in  the  vector  form  of  Eqn.  1  except  that  now  the 
various  terms  become: 


where  Rxx,  Rxy*  Ryx»  an<*  Ryy  are  the  same  as  before. 

In  addition  to  using  standard  approximate  factorization,  a  diagonally 
dominant  ADI  procedure  has  also  been  considered.  This  procedure  is 
slightly  more  stable  than  the  standard  ADI  procedure,  but  still  encounters 
instability  in  ranges  of  interest. 

Finally,  the  approximation  of  taking  pu  and  pv  as  known  quantities  and 
solving  the  temperature  equation  (Eqn.  3)  independently  and  In  combination 
with  the  radiation  equation  has  been  attempted.  This  procedure  works 
reasonably  well  and  provides  solutions  for  the  laser  absorption  problem. 

To  couple  this  equation  with  the  laser  equation,  we  express  the  laser 
equation  as, 

It  +  fi*6  "  '  Gk*  (5) 

where  8  =  in  Ij,A  represents  the  natural  logarithm  of  the  power  in  each 
ray  of  the  laser  beam,  c  represents  a  (fictitious)  speed  of  light,  and  k 
represents  the  absorptivity  of  the  medium.  When  Eqn.  5  and  3  are  coupled 
and  solved  simultaneously,  they  give  a  matrix  like  that  shown  in  Fig.  1. 
When  the  absorptivity  is  independent  of  temperature  and  pressure,  the  two 


equations  are  uncoupled  and  the  solution  can  be  found  by  solving  a 


i 

tridiagonal  matrix  followed  by  a  lower  triangular  matrix.  When  k  depends 
on  the  properties,  the  two  equations  are  coupled  and  generate  a  lower 
triangular  plus  tridiagonal  matrix  of  the  type  shown  in  Fig.  1.  Solution 
of  this  coupled  system  provides  convergence  that  is  as  rapid  with  laser 
heating  as  it  is  when  a  specified  heat  addition  rate  is  used.  That  is  to 
say,  when  the  energy  and  radiation  equations  are  solved  in  coupled  fashion, 

i 

the  laser  heat  addition  does  not  slow  convergence.  When  the  coupling 
between  these  two  equations  is  brought  in  iteratively,  laser  absorption 
causes  the  convergence  to  slow  down  dramatically  eventually  leading  to 
divergence.  These  results  show  that  coupling  between  the  laser  and  energy 
equations  is  a  desirable  aspect  of  any  numerical  solution.  Some 
representative  solutions  of  these  two  equations  are  described  in  the  next 
paragraph.  Algorithms  based  on  iterating  between  the  energy/radiation 
equation  system  and  the  continuity/momentum  equation  system  have  thus  far 
proved  unfruitful. 

Some  results  for  the  coupled  energy  and  laser  radiation  equations  are 
given  in  Figs.  2-9.  Figures  2  and  3  show  representatives  of  two  families 
of  grids  that  have  been  used.  Both  are  for  hollow  beam  laser  radiation  and 
the  geometric  shape  of  the  converging  beam  can  be  seen  by  the  rays  used  to 
represent  it.  The  fluid  dynamic  grid  has  been  treated  as  equally  spaced  as 
shown  on  Fig.  2,  or  highly  stretched  as  given  in  Fig.  3»  In  general, 
calculations  on  the  stretched  grid  proved  most  effective. 

Figures  H,  5  and  6  show  the  results  of  varying  inflow  velocity  for  a 
laser  power  of  720  W.  Figure  4  is  for  1.0  m/s,  Fig.  5  for  0.15  m/s,  and 
Fig.  6  for  0.01  m/s.  Figure  7  shows  a  close-up  view  of  temperature 
profiles  near  the  focal  point  for  the  case  shown  on  Fig.  6.  At  1.0  m/s 


(Fig.  4),  the  absorption  zone  is  pushed  back  almost  to  the  focal  region  of 
the  laser  beam.  As  the  velocity  is  reduced,  the  absorption  zone  remains 
almost  constant.  For  the  entire  range  of  velocities,  the  peak  temperature 
is  nearly  constant,  varying  from  17,800  K  at  the  high  speed  to  16,850  K  at 
the  low  speed. 

The  effect  of  changing  laser  power  at  an  inflow  velocity  of  0.01  m/s 
is  shown  on  Fig.  6  (720  W),  Fig.  8  (900  W) ,  and  Fig.  9  (1000  W).  For  this 
range  of  parameters,  changing  the  laser  power  has  little  effect  on  the 
shape  of  the  absorption  region,  and  the  peak  temperature  again  remains 
relatively  constant  although  it  does  decrease  slightly  as  power  is 
increased. 
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Power. 

Merkle,  C.  L.  and  Choi,  Y.-H.,  "A  Method  For  Computing  Compressible 
Flows  at  Very  Low  Mach  Numbers",  AIAA  Paper  86-0351,  AIAA  24th 
Aerospace  Sciences  Meeting,  Reno,  NV,  Jan.  1986.  To  appear  in  AIAA 
Journal,  May  1987. 

Merkle,  C.  L.  and  Choi,  Y.-H.,  "Computation  of  Low  Mach  Number  Flows 
With  Buoyancy",  Proceedings  of  the  10th  International  Conference  on 
Numerical  Methods  in  Fluid  Dynamics,  Beijing,  China,  June  23-27, 

1986,  pp.  169-173. 
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11.  Merkle,  C.  L.  and  Hosangadi,  A.,  "Computation  of  Low  Speed  Viscous 
Flows  with  Heat  Addition",  Sixth  IMACS  International  Symposium  on 
Computer  Methods  for  Partial  Differential  Equations,  Lehigh 
University,  June  23~26,  1987,  Bethlehem,  PA. 
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III.  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  RESEARCH  EFFORT 
Professional  Staff  - 

Charles  L.  Merkle,  Principal  Investigator,  Professor,  Mechanical 
Engineering 

Graduate  Students  - 

Michael  J.  Stanek,  Graduate  Assistant,  February  1981-August  1982. 
M.S.  Thesis  title,  "Analytical  Studies  of  the  Absorption  Mechanisms 
of  Equilibrium  Hydrogen."  Present  position:  Captain,  US  Air 
Force,  AFWAL/POTC,  WPAFB,  OH. 

Anil  Gulati,  Graduate  Assistant,  September  1981-August  1983. 

M.S.  Thesis  title,  "The  Absorption  of  Electromagnetic  Radiation  in 
an  Advanced  Propulsion  System."  Present  position:  Research 
Scientist,  GE  Corporate  Research  and  Development  Center, 
Schenectady,  NY. 

Gregory  A.  Molvik,  Graduate  Assistant,  September  1982-January  1985. 
M.S.  Thesis  title,  "A  Two-Dimensional  Analysis  of  Laser  Heat 
Addition  in  Converging  Nozzles".  Present  address:  CFD  Group, 
NASA/Ames,  Moffett  Field,  CA. 

Peter  Tsai,  Graduate  Assistant,  September  1983-August  1985. 

M.S.  Thesis  title,  "Stability  Characteristics  of  Laser-Supported 
Plasmas".  Present  position:  Ph.D.  Candidate  Penn  State  University. 

Yun-Ho  Choi,  Ph.D.  Candidate,  January  1985  to  present. 

Anticipated  Ph.D.  Thesis  title,  "Computation  of  Low  Mach  Number 
Flowfields  with  Strong  Heat  Addition". 

Eric  Shaw,  Graduate  Assistant,  September  1984  to  May  1986. 

M.S.  Thesis  title,  "Laser  Absorption  in  Flowing  Gas".  Present 
position,  Ph.D.  Candidate,  Pennsylvania  State  University. 

Ashvin  Hosangadi,  Ph.D.  Candidate,  September  1985-present. 
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IV.  INTERACTIONS/SPOKEN  PRESENTATIONS 


"Analysis  of  Laser-Supported  Combustion  Waves  in  Flowing  Media", 
AFOSR/AFRPL  Rocket  Propulsion  Research  Meeting,  Lancaster,  CA,  March 
26,  1981. 

"The  Potential  for  Using  Laser  Radiation  as  an  Energy  Source  for 
Propulsion",  Orbit  Raising  Propulsion  Workshop,  Orlando,  FL,  January 
16,  1982. 

"The  Use  of  Electromagnetic  Radiation  as  an  Energy  Source  for 
Propulsion",  Symposium  on  Advanced  Propulsion  Concepts  Using 
Time-Varying  Electromagnetic  Fields,  East  Lansing,  MI,  February  4, 
1982. 

"Analysis  of  Laser-Supported  Plasmas  in  Flowing  Media",  AFOSR/AFRPL 
Rocket  Propulsion  Research  Meeting,  Lancaster,  CA,  March  3,  1982. 

"The  Absorption  of  Electromagnetic  Radiation  in  an  Advanced  Propulsion 
System",  AIAA  Electric  Propulsion  Meeting,  New  Orleans,  LA,  November 
19,  1982. 

"Prediction  of  the  Flowfield  in  Laser  Propulsion  Devices",  AIAA  18th 
Thermophysics  Conference,  Montreal,  Canada,  May  1983. 

"The  Effect  of  Strong  Heat  Addition  on  the  Convergence  of  Implicit 
Schemes",  AIAA  6th  CFD  Conference,  Danvers,  MA,  June  1983. 

"A  Two-Dimensional  Analysis  of  Laser  Heat  Addition  in  Converging 
Nozzles",  AIAA  Aerospace  Sciences  Meeting,  Reno,  NV,  January  1984. 

"High  Power  Nd-Glass  Laser  Instrument  for  Advanced  Propulsion  and 
Diagnostics",  AFOSR/AFRPL  Rocket  Propulsion  Research  Meeting, 
Lancaster,  CA,  March  12-15,  1984. 

"Stability  of  Absorption  Phenomena  in  Laser-Thermal  Propulsion 
Gasdynamic  Interaction",  AIAA  17th  Fluid  Dynamics,  Plasmadynamics  & 
Lasers  Conference,  Snowmass,  CO,  June  25-27,  1984. 

"Analytical  Modeling  of  Strong  Radiation  Gasdynamic  Interaction", 
AFOSR/AFRPL  Chemical  Rocket  Research  Meeting,  Lancaster,  CA,  March  20, 
1985. 

"Numerical  Solution  of  Strong  Radiation  Gasdynamic  Interactions",  AIAA 
18th  Fluid  Dynamics,  Plasmadynamics  &  Lasers  Meeting,  Cincinnati,  OH, 
July  16-18,  1985. 

"An  Implicit  Time-Dependent  Scheme  for  Low  Speed  Flow  with  Buoyancy", 
22nd  Annual  Meeting  Society  of  Engineering  Science,  Inc.,  University 
Park,  PA,  Oct.  7-9,  1985. 

"Computation  of  Compressible  Flows  at  Very  Low  Mach  Numbers",  AIAA 
Aerospace  Sciences  Meeting,  Reno,  NV,  Jan.  6-9,  1986. 
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"Interaction  Between  Radiation  and  Gas  Dynamics",  invited  lecture, 
Chemical  Engineering  Department  Seminar  Series,  Michigan  State 
University,  East  Lansing,  MI,  February  10,  1986. 

"Computation  of  Low  Mach  Number  Flows  With  Buoyancy",  presented  at  the 
10th  International  Conference  on  Numerical  Methods  in  Fluid  Dynamics, 
Beijing,  China,  June  23-27,  1986. 

"Numerical  Solution  of  Strong  Radiation  Gas  Dynamic  Interactions", 
AFOSR/AFRPL  Chemical  Rocket  Research  Meeting,  Lancaster,  CA,  September 
1986. 

"Computation  of  Low  Speed  Viscous  Flows  with  Heat  Addition",  to  be 
presented  at  the  Sixth  IMACS  International  Symposium  on  Computer 
Methods  for  Partial  Differential  Equations,  Lehigh  University,  June 
23*26,  1987,  Bethlehem,  PA. 
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V.  INTERACTIONS/ADVISORY  FUNCTIONS 

Member  of  Workshop  Panel,  "Concepts  and  Experiments",  NASA/Michigan 
State  Symposium  on  Advanced  Propulsion  Concepts  Using  Time-Varying 
Electromagnetic  Fields,  February  1982. 

"Aerospace  Propulsion  at  Penn  State",  presentation  to  General  Robert 
T.  Marsh,  USAF,  The  Pennsylvania  State  University,  University  Park, 
PA,  May  5,  1982. 

"Modeling  of  Flowfields  that  Interact  with  Radiation  Fields", 
AFRPL/UDRI  Solar  Plasma  Propulsion  Workshop,  Bergamo  Center,  Dayton, 
OH,  Jan.  21-22,  1986. 


VI.  NEW  DISCOVERIES 


See  Publication  List. 
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Fig.  1  Tridiagonal-plus-lower  matrix  character  generated  when  laser 
equation  is  coupled  with  fluid  dynamics 


16 


1  :  M  l!  -t 

'-rtf-  i.M  ti‘ 

Mill  lit 

ihiin; 

...  !'  1  1 

.  ililil  Hi. 

|i  : . 

1- 

mill  hi. 

IMPl,  l; 

iliil 

mini  mi 

Imiiimi.  • 

'  ;,M  1  ‘ . 

,  .mw  mi 

1  Mi;  ,!'■■ 

.ii'min 

IIIIIIII  ill 

ll  Ililil 1 

:  1  I'll  ' 

-  i!  lillllll  III! 

I  iil:ini :tii!H 

I  :iii:!ii ; ;  i  r  :.i'i 

ll  Iiiiiiii  mm  'ililil 
ll  I'HIII!  !  Ill 

II  lillliil  hi  mi  i'll1 
|<  limilll.!:i,li  J.  -  i  ll 
ll  iilli;i; :  i:.1.:  it;  .1111 
|liiiiP;:ii;i  ill 
||n|:!!ll  .:;i  i:[ 

Ill.lPi  Wllilifi.  .'If 

ll  ri.il  r 
|lillllll;;!i!!li  ::MI'  it 

1 1 !  i  1 1 1 II I  (: :  i'  i  i  f I'  i  II 
||  IIIIIIII  ■l!!!il  . .  ■  I ' 
IMii-iiiii  liinn  i; 

ii  i  ■  ..ii  i  ■ 

|l  I  ;'i  |l: 
i.iiimminp 
pit  ■  ■ 

i  iii’ii  ..I  ;  . 

.1 

|li.;inii.i  i  i 

II  i, ;,ii  mi 

ll  Hill'll  !'!;ii!IC  . 
ICIIIIIill  IIH.1 
Iliililllll.i.I.II' 

I  ilium  iiiiiiii 

II  IIIIIIII  lllillL 

Iiihiiiii  mill'- 

limiiiiiiiiii'. 

Hill  IIP 
Ipiiium- 
liiiiiiiiiiit 
liiiiiiiiiit 


,:ii'l 
■i  "111. 

illini 
ill  IP 
Jiili!:iii 
"iinpr 
III  llll'll. 
mi  p,  li¬ 
nin';  in 

,:!lll  HI.I.IIl 
illlll  lllll'ili 
piiiiiiini; 
niiiiiniiiiiii 


,n  iiiiiiii  ml 

iiiiiiiiiiitiiil 

-mu  ilium  im| 

.lilliililllll 

iiiiiniiiiiiiii!i| 

iniii  iiiiiiii  ini 


'iiiiiiiiiiiiiiiiii; 

:  Pill  IIIIIIII  llll 
1  Illlll  IIIIIIII  III! 
Pl'llllllllllllllllll 
l.il  IIIIIIII  IIIIIIII  III! 
Il'il  IIIIIIII  IIIIIIII  llll 
:n , .  1 1 ! '.  I  i  IIIIIIII  llll 
i!  nil  lllllll  P 
Mill. :WIH. 
;ip-ii;, miii  i:;: 

•  „|:  "■'Iiiiiiii. 

.  I';,'  Ililil !'!, 

,  p;;i!'i:llil!ll! 

;  mi. illlll  ill: 

iiiiii  .mill  mi 
i  •.ipipiiiiiiii 
imi.iiiiiiliiii 

■iiiiii  ;iiiii!i  mil 
nun  iiiiiiii  mil 


7 Jr' 


V;VTw’'.*Vv. 


mmimr 


c 

• 

o 

CO 

o 

a 

B 

L. 

CO 

O 

c 

CO 

>» 

& 

•o 

CO 

TD 

0) 

3 

(0 

CO 

Cm 

i-H 

im 

c. 

o 

o 

TJ 

c 

a> 

0 

CO 

3 

4-> 

CO 

*o 

t. 

■H 

3 

bO 

to 

•o 

C 

a; 

o 

o 

o 

CO 

a 

B 

CO 

c0 

0) 

>» 

& 

£~ 

aj 

<V 

3 

00 

a* 

<0 

w 

» 

X3 

• 

c 

c 

CO 

o 

>1 

4-> 

u 

aj 

•P 

0) 

3 

B 

O 

o 

H 

<v 

CO 

o 

O 

CM 

00 

c*- 


lculation.  Stretched  grid  used  for  fluid  dynamics 


Temperature  contours  for  laser  power  of  720  W  and  inlet  velocity  of 
0.15  m/s.  Tmax  =  17,526  K;  AT  =  1000  K. 
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A.  Statement  of  Work 

Determine  experimentally  the  conditions  under  which  plasmas  can  be 
initiated  and  sustained  in  the  three  energy  addition  modes  (TMqi,  TEqi  , 
and  planar).  Initial  testing  will  be  with  nitrogen  and  helium. 

Nitrogen  will  be  used  to  simulate  the  molecular  nature  of  hydrogen 
while  helium  will  simulate  atomic  hydrogen.  Final  testing  will  be  with 
hydrogen.  Parameters  to  be  examined  include  gas  composition,  pressures 
and  flow  rates  and  microwave  power.  Quantify  system  heat  losses. 
Measure  spectroscopically  electron  and  ion  temperatures  and  densities 
in  microwave  generated  plasmas  in  the  three  energy  addition  modes.  Due 
to  high  gas  pressures,  local  thermodynamic  equilibrium  (LTE)  will  be 
initially  assumed  although  nonequilibrium  effects  will  be  examined. 

The  principal  measurement  will  be  ion  temperature  as  this  translates  to 
thrust.  The  effects  of  various  hot  and  cold  gas  mixing  schemes  on 
final  temperature  will  be  studied. 


B.  Status  of  Research  Effort 

Microwave  equipment  capable  of  generating  up  to  3  kW  CW  power  at  a 
frequency  of  2. 45  GHz  has  been  assembled  and  checked  out.  A  microwave 
resonant  cavity  wherein  both  filamentary  and  toroidal  plasmas  can  be 
produced  was  custom  manufactured  under  the  direction  of  Dr.  Jes  Asmussen  of 
Michigan  State  University  and  has  been  interfaced  with  the  microwave 
generator.  Additional  experimental  equipment  required  to  generate  planar 
propagating  plasmas  has  also  been  procured.  A  0.5  meter  Spex  scanning 
monochrometer  with  a  photomultiplier  tube  light  sensor  has  been  interfaced 
with  an  IBM  PC  to  allow  the  microcomputer  to  both  control  the  scanning  of 
the  monochrometer  and  the  recording  of  the  light  intensity  signal.  Quartz 
gas  containment  vessels  to  be  placed  inside  the  resonant  cavity  capable  of 
withstanding  up  to  3  atmospheres  of  internal  pressure  have  been  purchased 
and  the  associated  plumbing  to  flow  gases  over  a  range  of  flow  rates 
through  the  quartz  vessels  has  been  assembled.  The  production  of 
filamentary  microwave  heated  plasmas  in  nitrogen  gas  has  been  initiated. 
Figure  1  shown  the  three  microwave  absorption  modes  being  examined  in  this 
program.  Parameters  to  be  determined  are  which  combinations  of  gas 
composition,  pressure  and  flow  rate  will  permit  a  plasma  to  be  sustained  in 
the  resonant  cavity  for  a  given  microwave  power  input  and  standing  wave 
mode. 

A  numerical  model  of  the  one-dimensional  planar  propagating  microwave 
plasma  in  hydrogen  gas  had  previously  been  successfully  formulated.  This 
model  solves  the  system  of  governing  equations  consisting  of  the 


SCIENTIFIC  APPROACH 


FIGURE  1 


one-dimensional  steady  energy  equation 
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where  p  is  the  gas  density,  u  is  the  gas  velocity,  Cp  is  the  specific  heat, 
T  is  the  gas  temperature,  X  is  the  thermal  conductivity,  o  is  the 
electrical  conductivity  and  E  is  the  electric  field  vector;  and  Maxwell's 
equation  describing  the  propagation  of  the  microwave  energy 
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where  e  is  the  permit! vity  and  y  is  the  permeability.  The  specific  heat, 
thermal  conductivity  and  electrical  conductivity  are  functions  of 
temperature.  The  electromagnetic  energy  is  absorbed  by  the  plasma  as  the 
temperature  rises  with  some  microwave  energy  being  reflected  or  transmitted. 
Due  to  thermal  conduction  to  the  cold  gas  ahead  of  it,  the  plasma 
propagates  toward  the  microwave  energy  source  at  a  velocity  determined  by 
the  energy  balance  between  the  absorbed  microwave  power  and  the  heated  gas 
which  is  convected  away  downstream.  The  two  governing  equations  were 
numerically  integrated  using  a  fifth/sixth  order  variable  step  Runge-Kutta 
scheme.  An  iterative  method  was  used  to  determine  the  propagation  velocity 
eigenvalue,  pu,  similar  to  the  method  used  by  Kemp  and  Root  to  solve  for 
the  propagation  velocity  of  the  laser  heated  plasma. 

Since  then,  more  accurate  data  for  the  electrical  conductivity  of 
hydrogen  gas  at  elevated  temperatures  obtained  by  Yos  has  been  incorporated 
into  the  model  along  with  the  temperature  dependent  properties  of  helium 
gas.  Figure  2  plots  the  propagation  velocity  as  a  function  of  incident 
power  for  hydrogen  and  helium  at  one  atmosphere  pressure  showing  increased 
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plasma  propagation  velocity  with  increased  microwave  power.  It  can  also  be 
seen  that  for  the  same  microwave  power  level  the  plasma  propagates  in 
hydrogen  about  four  times  faster  than  in  helium.  Figure  3  plots  the 
maximum  gas  temperature  as  a  function  of  incident  power  and  shows  a  small 
increase  in  maximum  gas  temperature  as  the  microwave  power  is  increased. 
Figure  4  plots  the  percent  power  absorbed  by  the  propagating  plasma  as  a 
function  of  incident  power.  It  was  found  that  hydrogen  absorbs 
approximately  four  times  more  power  than  helium  and  that  for  both  gases  all 
of  the  incident  power  was  either  reflected  or  absorbed  with  no  power  being 
transmitted  through  the  plasma.  The  results  of  the  numerical  model  will  be 
compared  to  experimentally  measured  values  as  they  become  available. 
Calculations  are  planned  for  nitrogen  gas  and  for  all  three  gases  at  a 
pressure  of  ten  atmospheres.  Solutions  for  these  gases  are  desired  since 
initial  experimental  tests  will  use  both  nitrogen  and  helium.  Other 
planned  modifications  include  the  inclusion  of  radiation  losses  by  the  hot 
gas  and  the  examination  of  two-dimensional  effects.  Two-dimensional 
effects  are  important  because  microwave  radiation  transmitted  through  a 
waveguide  is  not  spatially  uniform  in  a  direction  normal  to  the  propagation 
direction.  Thus  previous  experimental  investigations  observed 
columnar-shaped  plasmas  instead  of  the  planar-shaped. 
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C.  Cumulative  Chronological  List  of  Written  Publications  in  Technical  Journals 


1.  Knecht,  J.  P.  and  Micci,  M.  M.,  "Analysis  of  a  Microwave-Heated  Planar 
Propagating  Hydrogen  Plasma."  Submitted  for  publication  in  AIAA 
Journal. 


2.  Durbin,  M.  R.  and  Micci,  M.  M.,  "Analysis  of  Propagating  Microwave 
Heated  Plasmas  in  Hydrogen,  Helium  and  Nitrogen,"  AIAA  Paper  87-1013. 
To  be  submitted  for  publication  in  AIAA  Journal. 
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D.  Professional  Personnel  Associated  with  Research  Effort 
Professional  Staff 

Michael  M.  Micci,  Assistant  Professor,  Aerospace  Engineering 
Graduate  Students 

James  P.  Knecht,  Graduate  Assistant,  January  1985-May  1986.  M.  S. 
Thesis  title,  "Numerical  Analysis  of  a  Microwave-Heated  Planar 
Propagating  Hydrogen  Plasma."  Present  position:  MIT  Lincoln  Labs, 
Lexington,  MA. 

William  Maul,  Graduate  Assistant,  June  1985-present.  Anticipated 
M.  S.  Thesis  title,  "Investigation  of  Resonant  Cavity  Microwave-Heated 
Plasmas." 

Michael  R.  Durbin,  Graduate  Assistant,  July  1986-present.  Anticipated 
Ph.D.  Thesis  title,  "Analytical  and  Experimental  Study  of  Propagating 
Microwave-Heated  Plasmas." 
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E.  Interactions 


(1)  Spoken  papers 


Coupling  Between  Gas  Dynamics  and  Microwave  Energy  Absorption. 
Presented  at  the  AFOSR/AFRPL  Rocket  Propulsion  Research  Meeting,  March 
21,  1985,  Lancaster,  CA. 

Prospects  for  Microwave  Heated  Propulsion.  Presented  at  Aerojet  Tech 
Systems,  Sacramento,  CA,  Oct.  8,  1985. 

Prospects  for  Microwave  Heated  Propulsion.  Presented  to  Project 
Forecast  II,  November  19,  1985,  Arlington,  VA. 

Coupling  Between  Gas  Dynamics  and  Microwave  Energy  Absorption. 
Presented  at  the  AFOSR/AFRPL  Rocket  Propulsion  Research  Meeting, 
September  11,  1986. 

Microwave  Electrothermal  Propulsion.  Presented  at  the  Air  Force 
Office  of  Scientific  Research,  Bolling  AFB,  DC,  September  18,  1986. 
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